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ABSTRACT 

A model of the ferromagnetic origin of magnetic fields of neutron stars is considered. In 
this model, the magnetic phase transition occurs inside the core of neutron stars soon 
after formation. However, owing to the high electrical conductivity the core magnetic 
field is initially fully screened. We study how this magnetic field emerges for an outside 
observer. After some time, the induced field that screens the ferromagnetic field decays 
enough to uncover a detectable fraction of the ferromagnetic field. We calculate the 
time-scale of decay of the screening field and study how it depends on the size of 
the ferromagnetic core. We find that the same fractional decay of the screening field 
occurs earlier for larger cores. We conjecture that weak fields of millisecond pulsars, 
B ~ 1 8 — 10 9 G, could be identified with ferromagnetic fields of unshielded fraction 
e ~ 10~ 4 — 10~ 3 resulting from the decay of screening fields by a factor 1 — e in ~ 10 8 
yr since their birth. 
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1 INTRODUCTION 

The physical origin of the magnetic field of neutron stars 
remains an open problem. Some researches adopt a working 
hypothesis that the field is inherited from the progenitor 
star. There also exist more specific models of generation of 
the magnetic field during the early cooling of a young neu- 
tron star. Generally, this kind of approach links the mag- 
netic field of neutron stars to macroscopic currents gener- 
ated either by the motion of the ionized material or by heat 
transport (Chanmugam 1992). 

For neutron stars, unlike for any other astrophysical ob- 
ject, there exists a different physical possibility, namely that 
the magnetic field is of microscopic origin and reflects spin- 
polarization properties of the neutron star matter itself. The 
ground state of dense baryonic matter is determined by nu- 
clear interactions that possess very strong spin components. 
It seems likely that these interactions lead to spin ordering 
in the neutron star matter making the ground state of the 
system a permanently magnetized phase, at least in some 
range of densities. 

Mechanisms of spin ordering in dense matter have been 
studied by a number of authors (Kutschera & Wojcik 1989, 
1990, 1993, 1996, 1997; Kutschera, Broniowski & Kotlorz 
1990a,b; Niembro et al. 1990; Marcos et al. 1991; Kutschera 
1994a,b). Kutschera & Wojcik (1989) have shown that neu- 
tron star matter with a low proton fraction of a few per 
cent is particularly susceptible to ferromagnetic spin order- 
ing. We have shown (Kutschera & Wojcik 1990, 1993, 1995; 



Kutschera 1994a,b) that protons in neutron star matter with 
such a small proton admixture behave as polarons that can 
localize at higher densities. Neutron star matter with local- 
ized protons displays a ferromagnetic instability (Kutschera 
& Wojcik 1989, 1996,1997) and the spin-ordered system de- 
velops a permanent magnetization, M ~ np/i e //, where 
np is the proton number density and fi e ff is the effective 
magnetic moment of the localized proton. Typical values of 
magnetization are in the range 10 13 — 10 14 G (Kutschera & 
Wojcik 1996). 

There exist also other model mechanisms that result 
in a spontaneous magnetization of hadronic matter inside 
neutron stars. The presence of the neutral pion condensate 
in the nucleon matter can lead to ordering of nucleon spins 
(Dautry & Nyman 1979). The same mechanism can operate 
in quark matter with chiral condensate (Kutschera et al. 
1990a, b). The magnetization can be higher in this case, up 
to ~ 10 15 G. 

Uncertainties in the physics of dense hadronic mat- 
ter prevent a firm conclusion concerning the nature of the 
ground state of baryon matter in neutron stars. However, 
ubiquity of spin interactions between hadrons with strong 
couplings makes the spin-polarized ground state of dense 
baryon matter a possibility that deserves serious considera- 
tion. 

Independently of the mechanism of polarization, one 
can study the phenomenological consequences of the exis- 
tence of spin-polarized baryon matter in neutron stars. If 
the ground state of dense matter is permanently magnetized 
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at baryon number densities exceeding some critical density 
n f , then there could exist a class of neutron stars possessing 
magnetic cores provided n/ < n™ ax , where n™ ax is the cen- 
tral baryon density of the maximum mass neutron star. Only 
sufficiently massive stars, with central densities n c > n/ 
would belong to this class. The critical density, n/, is esti- 
mated to exceed twice the nuclear saturation density. The 
magnetized core is a source of the magnetic field. However, 
the phase transition producing the magnetized core occurs 
in the inner part of the star, which has very high electrical 
conductivity. Thus, any magnetic field due to the sudden ap- 
pearance of the spontaneous polarization of baryonic matter 
becomes fully screened in such a medium. It takes some time 
for the magnetized core to become detectable by an outside 
observer. The aim of this paper is to investigate this prob- 
lem. 

The fact that any magnetic moment spontaneously cre- 
ated in the neutron star core is initially fully screened has 
important consequences for the spin structure of the mag- 
netic core, which can form a single magnetic domain. Such 
a domain is the most favourable configuration as far as the 
hadronic energy gain resulting from spin ordering is con- 
cerned. Because of screening, electromagnetic interactions 
do not force the single domain core to fragment into smaller 
magnetic domains. 

In the next section we describe the model of the magne- 
tized core. In Section 3 we study the screening of the ferro- 
magnetic field. Then in Section 4 the decay of the screening 
currents is considered. Section 5 contains astrophysical im- 
plications of the model. Some useful formulae are given in 
Appendix A. 

2 THE MAGNETIZED CORE INSIDE A 
NEUTRON STAR 

A single domain structure of the magnetic core is favoured 
by the dynamics of the magnetic phase transition. As the 
neutron star matter cools, the conditions for the sponta- 
neous spin ordering to occur are met for the first time at 
the centre of the star. This is because the energy gain re- 
sulting from spin ordering in the ferromagnetic phase in- 
creases with density (Kutschera & Wojcik 1990, 1993, 1995, 
1996). Thus, additionally, the critical temperature increases 
with density, and, for a roughly isothermal core, the first 
small bubble of the spin-ordered phase is nucleated at the 
centre of the neutron star. This first small domain grows as 
the temperature drops with nucleon spins in the next lay- 
ers of matter being polarized in the same direction. A single 
domain is the minimum hadronic energy configuration of 
the spin-polarized phase of neutron star matter. So far we 
have neglected electromagnetic interactions in the polarized 
phase. One can worry that, when magnetostatic energy is 
included, such a single domain may fragment into smaller 
randomly oriented domains. We show at the end of Section 

3 that this is not the case. 

To model the magnetized core we assume that the spon- 
taneously polarized matter forms a single domain with mag- 
netization in the z-direction: M = M(r)z. For this study the 
only relevant quantity is the magnetization as a function of 
baryon number density, M = M(tib). Calculations reported 
by Kutschera & Wojcik (1996) show that the magnetization 
depends only weakly on density, M w Mo ~ 10 13 - 10 14 G in 
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the whole range of relevant densities. In the case of polarized 
quark matter with chiral condensate, the magnetization is 
essentially density independent (Kutschera et al. 1990a,b). 
We can thus parametrize the magnetization distribution in- 
side neutron stars as 

M(r)=M f(r), (1) 

where the function f(r) accounts for the radial variation 
of magnetization. In the following we adopt a simple form: 
f(r) = 1 for r < r a and f(r) = for r > rv This function 
corresponds to a uniformly magnetized sphere of radius r a 
with a surface layer of thickness h = — r a where the 
magnetization drops linearly to zero at r = n with f(r) = 
(r — rt)/(r a — 77,) (Fig.l). This simple function is flexible 
enough to account for the main features of the magnetization 
distribution in neutron stars. 

The magnetic moment of the core is 

d cora = / M{r)d 3 r, (2) 
Jo 

with M(r) given by equation (1). The size of the core is 
determined by the critical density for spontaneous polariza- 
tion, n/ . Thus the nonzero magnetization, M(r) / 0, ex- 
ists only in the inner core of stars of density n_g(r) > n/. 
The core radius r b corresponds to the critical density nf. 
ns(n>) = nf. 

The contribution of the magnetic moment (2) to the 
magnetic field at the pole is 



where Rns is the radius of the neutron star. This field is 
Bp ~ 10 12 — 10 la G for the above typical values of the 
magnetization. 

3 SCREENING OF THE FERROMAGNETIC 
FIELD 

Calculations in various models suggest the energy per 
baryon in the polarized phase could be below that for the 
normal phase by at least ~ lMeV. The phase transition 
from normal matter, with no spin ordering, to magnetized 
matter is thus expected to occur very soon after the forma- 
tion of the neutron star. The phase transition is completed 
quickly. For any practical purposes one can assume that the 
ferromagnetic field is switched on instantaneously. This is 
because the time-scale for magnetic field diffusion, deter- 
mined by the electrical conductivity, is much longer than 
the duration of the phase transition. 

The neutron star matter outside the magnetic core is 
a medium of very high electrical conductivity a. Realistic 
calculations show that conductivities corresponding to the 
neutron star crust are, typically, <j crust ~ 10 23 s _1 (Chan- 
mugam 1992). These values are lower than those correspond- 
ing to the core, a core ~ 10 29 s _1 (Chanmugam 1992). It is 
rather obvious that sudden switching-on of the magnetic 
field of the magnetized core will result in the induction of 
the screening field, which will fully shield the ferromagnetic 
field. This happens because for such a high conductivity the 
magnetic flux through any loop is conserved on the time- 
scale of the phase transition. Flux conservation then requires 
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Figure 1. The magnetization profile function f(x) and the func- 
tion g(x) for z b = 0.9. Solid and dashed curves correspond, re- 
spectively, to z a = 0.001 and z a = 0.8. 

that dJi/dt — 0. We assume that there is no magnetic field 
before the phase transition occurs, B = 0. Thus, at the in- 
stant t = the switched-on ferromagnetic field, B/ er , and 
the induced field, Bj rat j(t), cancel one another exactly, 



B fer + B ind (0) = 0. 



(4) 



The ferromagnetic field of the core is easily calculated. 
The vector potential has only the (^-component 



A$ — sinO, 



(5) 



where the function R(r) reads 

R(r)=d(r)/r. (6) 
Here d(r) is the magnetic moment inside the sphere of radius 



d(r) = I M(r')d 3 r'. 



B 



Components of the magnetic field are 
- ^cosO, 



(7) 



(8) 



sinO. 



(9) 



Bfer = 1 dR(r) 
6 r dr 

The flux conservation condition, equation (4) , gives the 
components of the induced field that screens the ferromag- 
netic field, 



B l r nd (0) = -B. 



fer 
r 5 



fer 



(10) 



(11) 



Be" (0) = -B i 

These formulae show that the induced field components, 
B r nd and B' g " d , satisfy equations (8) and (9) with the func- 
tion Rind{r,t), which, however, depends also on time. At 
t = the condition (4) implies that 



Rind(r,0) = -R(r). 



(12) 



One can also obtain the space structure of the current 
sustaining the induced field. Generally, this current changes 
in time, Ji„ d = 3i n d(r,0,t). From the formula V x B;„d = 
47r/cJ in£ j it follows that 



4w ji n d _ sin9 d Rj nd _ 2R in d . 

At t = using equation (4) , we find 
3 lnd (r,6,0) = -cV x M. 



(13) 



(14) 



The above discussion allows us now to address the ques- 
tion of the stability of a large single domain against fragmen- 
tation into smaller randomly oriented domains. In the case of 



terrestrial ferromagnets, large domains are not stable. Their 
magnetostatic energy can be lowered by producing a number 
of smaller domains with apparently randomly oriented mag- 
netic moments. Let us stress that reduction of the energy is 
due to dipole magnetic interactions between these magnetic 
moments. The magnetic core in neutron stars forming a sin- 
gle domain with screened magnetic field cannot lower its 
magnetostatic energy by fragmenting into smaller domains. 
This is because the magnetic moment of any new domain 
would also be fully screened and thus, unlike in terrestrial 
ferromagnets, there would be no magnetic dipole interac- 
tions between neighbouring domains. The core is thus stable 
against fragmentation into small domains. 



4 DECAY OF SCREENING FIELDS 

The induced current, J in d, will suffer ohmic decay as the 
electrical conductivity, a, though very high, is finite. The 
nonzero net field, B = B/ er + B ind (t) / 0, will eventually 
emerge. Let us note that ohmic decay is the only relevant 
mechanism of magnetic field decay as long as the net mag- 
netic field is low. Ambipolar diffusion and Hall drift, which 
play a crucial role in the decay of strong fields (Goldreich & 
Reisenegger 1992), can be safely neglected. 

To calculate the time behaviour of the screening field 
we apply a standard analysis of decay modes (Wendell, Van 
Horn & Sargent 1987). To avoid unnecessary complications 
in our exploratory study we neglect spatial and temporal 
dependence of the conductivity and assume that a = const. 
The time dependence of the induced field can be found from 
the expansion of the function Ri nd (r,t) into eigenfunctions 



X n (x) — V2nirxji(nnx), 



(15) 



where ji is the Bessel function and x — rj R-* is the normal- 
ized radial variable with a suitably chosen radius R*. The 
expansion reads 



Ri„d(x,t) = ^ C„X n (x) exp (-t/r n ), 



(16) 



where r n = 4ir R^a / '(am) 2 is the decay time of the nth 
mode. 

The expansion coefficients, C n , are obtained using the 
function Ri n d(x,0) — —R(x) as the initial condition: 



Cn = / Rind(x,0)X„(x)dx. 



(17) 



With our choice of the magnetization profile, f(r), the 
function R(x) reads 



R(x) = -TtM Rlg(x), 



(18) 



where the function g(x) is determined entirely by the magne- 
tization profile function f(r). The function g(x), calculated 
by integrating the magnetization in equation (7) , is given in 
Appendix A. 

The form (18) of the function R(x) indicates that the 
time behaviour of the screening field is sensitive to the mag- 
netization profile. To study this behaviour we have calcu- 
lated the coefficients C„ of the expansion, equation (17), 
which are also given in Appendix A. The unshielded frac- 
tion, e(x,t), of the ferromagnetic field emerging after time t 
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e(x,t) = 



|B/ er + Bj nc ((t)| 



B 



= 1 - 



/er| 



Rind(x,t) 
Rind{x,0) 



(19) 



We are interested mostly in the value of e(x, t) at the surface 
of the star, x = 1, which we denote as e(i) = e(l,i). 

The formula (19) shows that e(t) does not depend on 
the magnetization Mo- This means that the relative rate of 
emergence of the ferromagnetic field depends only on the 
geometry of the magnetized core. For the case of uniform 
conductivity we consider here this behaviour is even more 
universal. We can express the time behaviour of the emerg- 
ing magnetic field in terms of the dimensionless variable 
y = t/ri, where 



Tl = 



4iCcr 



(20) 



is the decay time of the longest-living (fundamental) decay 
mode. The formula (19) becomes 



<y) = i 



J2 n C n X n (l)exp(-yn 2 ) 



(21) 



The main results of our analysis are presented in Fig. 2 
where we show e(y) for indicated values of the parameters 
z a = r a /R* and Zb = rb/R* governing the spatial distri- 
bution of magnetization in our simple model. As one can 
notice, the shape of the curves is determined essentially by 
the value of Zb with much smaller influence of z a - 

The general tendency is that the higher is Zb the higher 
is e. The curves for the same Zb and various z a are very 
close to each other. We see that the time dependence of the 
unshielded fraction e is most sensitive to the total radial 
extension of the ferromagnetic core. The radial distribution 
of magnetization, which is our model is controlled by z a , 
has a much smaller influence on e. This means that the fer- 
romagnetic field for cores of larger extension becomes un- 
covered earlier than for smaller ones. The curves shown in 
Fig. 2 prove that the unshielded fraction e is very sensitive 
to Zb and is much less sensitive to z a . For each value of z b 
we show curves corresponding to two extreme values of z a : 
z a — 0.001 (dashed curves) and z a = 0.9zb (solid curves). 

In our discussion above we have made a number of sim- 
plifying assumptions regarding the neutron star structure. 
The conductivity is assumed to be uniform and constant 
in time, which is equivalent to constant temperature. Also, 
Euclidean geometry is used. These simplifications are not 
expected to affect our qualitative conclusions as far as the 
decay time-scale of the screening field is concerned. 



5 CONCLUSIONS AND IMPLICATIONS 

The main conclusion from the above analysis is that the 
emergence time of the ferromagnetic field, at low values of 
the unshielded fraction e, is very sensitive to the size of the 
polarized core. For large cores, comprising ~ 90 per cent of 
the radius, a fraction e ~ 10~ 4 of the core field is visible after 
t ~ 10~ 4 ti. For small cores, Zb ~ 0.4, the same fractional 
field is visible after t ~ 10 _1 Ti, i.e. a factor of ~ 10 3 later. 

To explore the astrophysical consequences of this model 
we must specify the relevant values of the decay time n , or 
equivalently the electrical conductivity a. For magnetized 
cores in neutron stars, the relevant conductivity is that of 



Figure 2. The unshielded fraction e as a function of the variable 
y = t/Vi for indicated values of the magnetized core parameters 
z a and Zb- Curves labelled 1,2 and 3 correspond to Zb ~ 0.99, 
Zb ~ 0.9 and Zb ~ 0.4, respectively. 



the liquid core matter, a c 



lO 2 



The corresponding 



decay time is n ~ 10 12 yr. One should notice that the choice 
of the o"core implies that the presence of the neutron star 
crust becomes somewhat irrelevant for our considerations 
here. This is because the crust conductivity is lower by a 
factor ~ 1CF 6 . An important consequence is that the value of 
the radius R* we use in the definition of the variable x should 
be identified with the radius of the liquid core, R, = Ri c , 
rather than the total neutron star radius (which includes 
the crust). Correspondingly, z a and Zb are the fractions of 
the liquid core radius, z a = r a /Ri c and z b = r b /Ri c . With 
this observation the values z b > 0.9 do not seem unrealistic, 
especially for soft equations of state. 

Applying our model to neutron stars we find that the 
magnetic field of the magnetized core formed soon after 
the birth of the star emerges in ~ 10 s yr at a level of 
~ 10~ 4 Bp er ~ 10 8 G provided the ferromagnetic core is large 
enough. Both the time-scale of 10 8 yr and the magnetic field 
of 10 8 G are typical for millisecond pulsars. One can thus 
conjecture that the magnetic fields of millisecond pulsars 
are due to spin-polarized matter inside neutron stars. 

The ferromagnetic origin of the magnetic fields of mil- 
lisecond pulsars could explain the discrepancy between the 
birth rate of low mass X-ray binaries (LMXB), which are 
supposed to be the progenitors of millisecond pulsars, and 
the number of millisecond pulsars in the Galaxy (Bhat- 
tacharya 1995). In this case, single neutron stars born with 
short rotation periods and low magnetic fields (much less 
than 10 7 G) could become millisecond pulsars in ~ 10 8 yr. 
This would also help to understand the existence of single 
millisecond pulsars without invoking the companion evapo- 
ration scenario. Also, this model would shed some new light 
on the old question of the decay time of the inherited mag- 
netic field of normal radio pulsars. Currently, a popular view 
is that only magnetic fields of neutron stars that accreted 
matter in binary evolution decay significantly. In the fer- 
romagnetic model, the inherited magnetic field of an iso- 
lated neutron star could decay in ~ 10 7 yr, as was found in 
early studies of the subject (Ostriker & Gunn 1969). There 
would be no need for direct connection between the amount 
of matter accreted by the neutron star in LMXBs and the 
amount of decay of its magnetic field. Recent analysis by Wi- 
jers (1997) provides some evidence against field decay being 
proportional to the mass accreted by a neutron star during 
the LMXB evolution, supporting our conjecture. 

Let us stress finally that looking for the presence (or ab- 
sence) of a ferromagnetic component of the magnetic field 
of neutron stars is of great importance for the physics of 
hadronic matter. Physical conditions prevailing in neutron 
stars (high density, low temperature and (3 — equilibrium) 
are not accessible to any laboratory experiment. Magnetic 
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field could serve as a direct probe of dense hadronic matter in 
the neutron star interior. Evidence in favour of the presence 
of spin-polarized matter inside neutron stars would intro- 
duce new qualitative features of the ground state of dense 
hadronic matter. An opposite conclusion would mainly con- 
strain the critical density n/. The fact that a magnetic field 
resulting from the ferromagnetic core grows with time, while 
the inherited field decays in time, could potentially allow one 
to distinguish the two components. 
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and, 



C = 



-0.75 



Z b - Z a 

and, finally, 

g(x) = —, z b <x<l, 



(A5) 



(A6) 



where 



1 ,3 



D = z a -\ (-(z - z b ) + z b (z b - z a )) 



Z b - Z a 4 

The expansion coefficients C„ in equation (19) read 



(A7) 



C„ = ~-7vM R t c n , 



(A8) 



where 

c n = 2[a n (z a ) + A(b n (z a ) - b n (z b )) + B(a n (z b ) - a n (z a )) 



C_ 

(wr) 

Here 



+ t „^ A {d n {z b ) - d„(z a )) + Db„(z b )]. 



(A9) 



a n (z) = 



-3mvzcos(n-Kz) + 3sin(nnz) 



(ni) 3 



(n-Kz) 2 sin(n-Kz) 
(nn) 3 

, . , sin(n-Kz) 
b n (z) = — f '-, 



(A10) 
(AH) 



and 

d„(z) = 8cos(nnz) — A{rm z) 2 cos(n-K z) + 8nnzsin(niTz) 



— (nnz) sin(nnz). 



(A12) 
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APPENDIX A: 

For our choice of the magnetization profile, f(r), the func- 
tion g(x) is a smooth function, which is given analytically 
below: 



g(x) = X , X < z a , 



(Al) 



A 2 ■', 

g(x) — — + Bx + Cx , Za < X < Zb, 

where 

3 , 0.754 - z b zl 

A — Z a + ■ 



Zb — Z a 



B = 



Zb 

Zb — Z a ' 



(A2) 

(A3) 
(A4) 



